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ABSTRACT 

Aims. In a recent study, we derived individual distances for 109 pre-main sequence stars that define the Lupus kinematic association of 
young stars. Here, we use these new distances to derive the masses and ages of Lupus T Tauri stars with the aim of better constraining 
the lifetime of their circumstellar disks. 

Methods. Using the photometric and spectroscopic information available in the literature, we computed the photospheric luminosity 
of 92 T Tauri stars in the Lupus association. Then, we estimated their masses and ages from theoretical evolutionary models. Based 
on Monte Carlo simulations and statistical tests, we compare the mass and age distribution of the classical T Tauri stars (CTTS) and 
weak-line T Tauri (WTTS) in our sample. 

Results. We show that the CTTSs are on average younger than the WTTSs and that the probability that both T Tauri subclasses are 
drawn from the same mass and age parental distribution is very low. Our results favor the scenario proposed earlier for the Taurus- 
Auriga association, where the CTTSs evolve into WTTSs when their disks are fully accreted by the star. Based on an empirical disk 
model, we find that the average disk lifetime for the T Tauri stars in the Lupus association is = 3 x 10® (M./Mg)® yr. 
Conclusions. We find evidence that the average lifetime of the circumstellar disks in the Lupus association is shorter than in the 
Taurus-Auriga association and discuss the implications of this result. 
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1. Introduction 


T Tauri stars (TTSs) are late-type low-mass (M < 2Mq) pre- 
main sequence stars that are usually associated with molecular 
clouds. They tend to be found in young associations (Joy 1945; 
Herbig 1962; Bertout 1989). They are classified into one of two 
types: (i) the classical T Tauri stars (CTTSs) that show strong 
emission lines indicative of accretion from a circumstellar disk, 
and (ii) the weak-line T Tauri stars (WTTSs) that exhibit en¬ 
hanced magnetic activity and show no signs of accretion. To ex¬ 
plain the existence of both subgroups, it has been hypothesized 
for a long time (Walter et al. 1988; Bertout 1989) that the CTTSs 
become WTTSs when their disks dissipate. Although this is not 
a problem for the older population (f > 10 Myr) of WTTSs, 
this picture implies that the timescale of disk dissipation should 
be a few Myr to justify the CTTSs and WTTSs coexisting (see, 
e.g., Fig. 15 of Kenyon & Hartmann 1995) in the Hertzsprung- 
Russell diagram (HRD). Indeed, previous studies (Haisch et al. 
2001; Takagi et al. 2014) indicate that the disk lifetime in differ¬ 
ent clusters is typically in the range 3-6 Myr. However, Fedele 
et al. (2010) demonstrate in a more detailed study that mass ac¬ 
cretion and dust dissipation in protoplanetary disks occur on a 
typical timescale of 2.3-3.0 Myr (see Fig. 4 of their paper), which 
implies a rapid disk evolution. 


Some progress in determining the lifetime of circumstellar 
disks of TTSs was achieved by Bertout et al. (2007), who used 
the kinematic parallaxes derived from a kinematic study of the 
Taurus-Auriga T association (Bertout & Genova 2006) and cal¬ 
culated the stellar luminosities, masses, and ages to readdress 
the question of the relationship between CTTSs and WTTSs. 


As a result, they obtained, for the first time, the lifetime of 
a circumstellar disk in terms of the mass of the parent star, 
Trf = 4 X 10® (M/Mq)®-^® yr, and demonstrated that the observed 
mass and age distributions of the TTSs in the Taurus-Auriga as¬ 
sociation can be explained by assuming that the CTTSs evolve 
into WTTSs when their disks are fully accreted by the star. 

While supporting the current scenario of disk evolution for 
the Taurus-Auriga T association, this first result needs to be con¬ 
firmed and compared to other star-forming regions (SFRs). One 
reason for extending this study to other SFRs is that the origin of 
the initial mass function remains unclear, but is expected to de¬ 
pend on the intrinsic properties of each SFR, such as clustering, 
mass segregation, and binarity (Bonnell et al. 2007). Moreover, 
the evolution of circumstellar (protoplanetary) disks depends on 
environmental criteria (Rosotti et al. 2014), and a comparison 
with other SFRs could provide important clues for constraining 
the timescale of planet formation in different environments. For 
example, Marinas et al. (2013) conclude that the dependence of 
disk fraction on stellar mass in Taurus and the NGC 2264 open 
cluster is different (see Fig. 10 of their paper), and they argue 
that the lower disk fractions in NGC 2264 might be due to the 
much higher stellar density than in Taurus. 

In the following, we perform a study similar to the one by 
Bertout & Genova (2006) for the TTSs in the Lupus SFR. The 
Lupus complex of molecular clouds is located in an environment 
surrounded by the more massive stars from the nearby Scorpius- 
Centaurus OB association (de Zeeuw et al. 1999). This makes 
the Lupus SFR an interesting target for comparative studies with 
Taurus-Auriga, because of the different environmental and dy¬ 
namical effects (e.g., solar winds and UV radiation) caused by 
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the existence of hot stars in this region. In a recent work, Galli 
et al. (2013, hereafter Paper I) investigated the kinematics of the 
Lupus SFR and performed a membership analysis. They iden¬ 
tified 109 pre-main sequence stars that define the comoving as¬ 
sociation of the Lupus SFR (hereafter, the Lupus association) 
based on a new convergent point search method (Galli et al. 
2012) and a modified version of the k-NN method. Then, they 
used stellar proper motions and radial velocities to calculate in¬ 
dividual distances for all group members. In this paper, we use 
these newly derived kinematic distances to accurately determine 
the masses and ages of the TTSs in the Lupus association and 
constrain the lifetime of their circumstellar disks. 


2. Physical properties of the comoving stars 

The sample of stars used in this study is based on the 109 co¬ 
moving stars of the Lupus association identified in Paper 1. One 
star, SSTc2dJ154013.7-340142, which was included in our sam¬ 
ple as a YSO candidate member of the Lupus clouds (see Com- 
eron et al. 2009), is excluded from the upcoming analysis, be¬ 
cause it turned out to be a background giant star (see Comeron 
et al. 2013). Our sample consists of 38 CTTSs (or CTTSs sys¬ 
tems), 67 WTTSs (or WTTS systems), and 3 Herbig Ae/Be stars 
(HAeBe). These numbers correct those given in Sect. 7.1 of Pa¬ 
per I, where some objects were either misidentified or have been 
reclassified in the meantime. The results and conclusions given 
in Paper I are obviously not affected by this minor correction. 
We restricted our analysis in subsequent sections to the TTSs in 
the Lupus association (i.e., CTTSs and WTTSs) and excluded 
the HAeBes in our sample that appear in small numbers and are 
not statistically significant for this study. 

The TTS subclass (CTTS or WTTS) from the stars in our 
sample were taken from the literature (see Table 1 for refer¬ 
ences), and no attempt has been made to reclassify them. How¬ 
ever, we compiled the Ha equivalent widths (EW) of our targets 
and calculated the near-infrared excess emission that is indica¬ 
tive of optically thick accretion disks to investigate whether the 
given TTS subclass was consistent with other diagnostics. We 
measure the excess emission by 

( 1 ) 

where (H-K) is the dereddened color, and (H-K)q the intrinsic 
color of the star (see Hughes et al. 1994; Wichmann et al. 1997a). 

We used the standard limit of lOA (see, e.g., Appenzeller & 
Mundt 1989) to distinguish between CTTSs (EW(Ha) > lOA) 
and WTTSs (EW(Ha) < lOA), together with the limit of 
A(H - K) > 0.09 found by Edwards et al. (1993) for a sam¬ 
ple of TTSs that show spectroscopic evidence of accretion. Eig- 
ure 1 (upper panel) illustrates our analysis and clearly separates 
the accreting TTSs (i.e., CTTSs) from the non-accreting TTSs 
(i.e., WTTSs). Alternatively, White & Basri (2003) use the Ela 
EW with a different cutoff that depends on the spectral type 
to distinguish between CTTSs and WTTSs. They propose that 
a TTS is indeed a CTTS if EW(Ha) > 3A for K0-K5 stars, 
EW(Ha) > lOA for K7-M2.5 stars, EW(Ha) > 20A for M3- 
M5.5 stars, and EW(Ha) > 40A for M6-M7.5 stars. We ap¬ 
plied this classification scheme to the stars with spectral type 
later than KO in our sample as shown in Eigure 1 (lower panel). 
We modified the first spectral type interval slightly to include 
two K6 stars in our sample and provide a TTS subclass to all 
stars. Our analysis presented in Eig. 1 not only confirms the 
TTS subclass found in several sources in the literature, but also 


gives us more confidence when discussing the different proper¬ 
ties of these two subgroups in the rest of this paper. Only one 
star (Sz 76) would be misidentified by this alternative classifi¬ 
cation (see Eig. 1), because its EW(Ha) - 10.3A (see Hughes 
et al. 1994) is marginally above the cutoff limit for its spectral 
type. However, Sz 76 has been classified recently as a WTTS by 
Wahhaj et al. (2010) and a Class III source by Lopez Marti et al. 
(2011). Moreover, its excess emission A(H-K) is clearly beyond 
the cutoff limit of 0.09 between accreting and non-accreting 
TTSs (see Edwards et al. 1993). We therefore retain Sz 76 as 
a WTTSs in our study. 

The first step toward accurately computing the photospheric 
luminosity of stars is to derive their visual extinction. We com¬ 
puted the stellar extinction Ay from the color excess in (V - Ic), 
(V - Rc), and (J - K) for the WTTSs in our sample, follow¬ 
ing the procedure described in Sect. 4.2 of Kenyon & Hartmann 
(1995). We set the extinction to zero when the derived Ay es¬ 
timates yield a non-physical negative value or a non-significant 
result (i.e.. Ay < (TAy). We took the mean of all non-zero Ay 
estimates as the final result with the standard deviation as the 
formal uncertainty. We used a total-to-selective extinction ratio 
Ry equal to 3.1 (Savage & Mathis 1979) and adopted the redden¬ 
ing relations for the various colors from Bessell & Brett (1988). 
The extinction for the CTTSs in our sample is estimated from the 
(Rc - Ic) color excess with Ay = 4.76 ■ [(Rc - Ic) - (Rc - 7c)o] 
as described in Sect. 2 of Cieza et al. (2005). 

To gain more confidence in our results, we decided to com¬ 
pute our Ay estimates again by using the same dereddening 
method for CTTSs and WTTSs. In this alternative approach we 
used only the optical color excess in (V - Ic), (V - Rc), and 
(Rc - Ic) for all stars in the sample (i.e., CTTSs and WTTSs). 
We concluded that both strategies yield compatible results and 
that our choice of using different methods to deredden CTTSs 
and WTTSs (as described above) has no significant impact later 
on the determination of stellar masses and ages. 

We used the individual parallaxes given in Tables 6 and 7 of 
Paper I to derive the photospheric luminosity of the Lupus co¬ 
moving stars. The stellar luminosities for the CTTSs in our sam¬ 
ple were calculated from the Ic band, because the standard pro¬ 
cedure of using the J flux systematically overestimates their lu¬ 
minosities (see discussion in Cieza et al. 2005). Eor the WTTSs, 
we derived the stellar luminosities from the V flux when the Ic 
flux was not available. We verified that both procedures yield 
equivalent results within 5% when both measurements are avail¬ 
able. We used the intrinsic colors, bolometric corrections, and 
temperatures for each spectral type given by Kenyon & Hart¬ 
mann (1995). The uncertainties on the derived luminosities take 
the various sources of photometric errors into account, but they 
are strongly dominated by the error budget of parallaxes. No op¬ 
tical photometry is available for two of our WTTSs (RXJ1546.7- 
3459 and RXJ1612.3-4012). We thus estimate their luminosity 
from the J flux given by the 2MASS catalog (Cutri et al. 2003). 
Another star, Sz 102, was removed from the sample because 
it is believed to be associated with bipolar jets and an edge-on 
disk that makes our Av/ estimate very imprecise (Krautter 1986; 
Hughes et al. 1994). Moreover, its spectral type is not well- 
defined, and our analysis underestimates the luminosity, putting 
it below the zero-age main sequence (ZAMS). 

To better compare our results with those derived for the 
Taurus-Auriga association, we used the same completeness lim¬ 
its as Bertout et al. (2007) by removing the stars (or stellar 
components) that have derived luminosities lower than 0.15Lo 
and spectral types later than M4. These values correspond to 
y ^ 16 mag and T^ff ^ 3300 K. The final sample of stars 
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Fig. 1. Plots of the Htt EWs against the infrared excess emission 
A(H - K) (upper panel) and spectral type (lower panel) to distinguish 
between CTTSs and WTTSs. 


with known photometry and spectral type that fulfill the above- 
mentioned criteria contains 30 CTTSs and 62 WTTSs. Physical 
properties of the Lupus stars discussed in this section are given 
in Table 1. Eighteen stars in our sample were identified as bi¬ 
nary systems in the literature. The photometry of these stars was 
corrected for binarity with the information (mainly the flux ratio 
of the components) provided by Reipurth & Zinnecker (1993), 
Brandner et al. (1996), Ghez et al. (1997), Guenther et al. (2001), 
Tori'es et al. (2006), and Romero et al. (2012). 


3. Revisiting the CTTSs and WTTSs in Lupus 

In the following we use the stellar luminosities and effective tem¬ 
peratures of the stars to place them in the HRD and discuss the 
age and mass distributions of the Lupus association. We derive 
the mass and age of the stars by interpolating between the evo¬ 
lutionary tracks and isochrones of the pre-main sequence stars 
models computed by Siess et al. (2000). Among the various 
models that exist in the literature, we chose the Siess et al. (2000) 
grid of pre-main sequence stars to allow for a direct comparison 
with the results in the Taurus-Auriga association (see Sect. 5) 
obtained by Bertout et al. (2007). The uncertainties affecting 
our mass and age estimates are estimated by defining the upper 
and lower limits for the stellar luminosities and effective tem¬ 
peratures within their error bars, and interpolating between the 
tracks and isochrones. We assumed an error of + 100 K on the ef¬ 
fective temperatures that roughly corresponds to an error of one 
subtype of spectral type for the late-type stars in our sample. Our 
mass and age results are presented in Table 1. 

The HRD of the Lupus association illustrated in Fig. 2 shows 
that a few WTTSs and CTTSs share the same region of the HRD. 
As already discussed in Paper I (see, e.g.. Fig. 15 of that pa¬ 
per), the WTTSs in our sample comprise an on-cloud population, 
which is located in the direction of the main star-forming clouds 
(Lupus 1-4) of the Lupus SFR, and a more dispersed off-cloud 
population surrounding the clouds. The question arises whether 
the two populations of WTTSs and the CTTSs of the Lupus asso¬ 
ciation exhibit the same mass and age distributions. Table 2 lists 
the average mass and age of each population for comparison and 
reveals the different properties of the CTTSs and WTTSs in our 
sample. 

Table 2. Mean masses and ages of the TTSs in the Lupus association 
derived from the Siess et al. (2000) models. 


Sample 

Stars 

Mean 

M/Mo 

Median 

M/Mo 

Mean 

log/ 

(t in Myr) 

Median 

log/ 

(t in Myr) 

CTTS 

30 

0.54 ± 0.04 

0.48 

6.32 ± 0.24 

6.28 

WTTS (on-cloud) 

33 

0.69 ± 0.07 

0.59 

7.05 ± 0.67 

6.55 

WTTS (off-cloud) 

29 

1.09 ±0.07 

1.10 

7.19 ±0.28 

7.11 

WTTS (all stars) 

62 

0.88 ± 0.05 

0.81 

7.12 ±0.48 

6.87 


We computed our mass and age estimates again by using 
the D’Antona & Mazzitelli (1997) tracks and the Baraffe et al. 
(1998) models to demonstrate that the different age and mass 
distributions are not an artifact from our choice of evolutionary 
model (see Table 3). Although evolutionary models computed 
by different groups yield different mass and age estimates for a 
given star, we note from Tables 2 and 3 that the results obtained 
with the various models employed in this work are consistent be¬ 
tween themselves. This allows us to confirm that (i) the CTTSs 
are on average younger than the WTTSs, and (ii) the off-cloud 
stars in the Lupus association are on average older than the on- 
cloud stars. 

We performed a Kolmogorov-Smirnov (KS) and a Wilcoxon 
rank sum test to investigate the null hypothesis that the age and 
mass distributions of the CTTSs and WTTSs in the Lupus as¬ 
sociation derive from the same distribution. The results of this 
analysis are given in Tables 4 and 5. Adopting the significance 
level of a = 0.05 (see, e.g., Feigelson & Jogesh Babu 2012), we 
conclude that the null hypothesis can be rejected in each case 
regarding stellar ages (see Table 4), which means that the age 
distributions of the CTTSs and WTTSs in Lupus are not identi¬ 
cal. On the other hand, our results with the statistical tests ap¬ 
plied to the mass distribution of CTTSs and on-cloud WTTSs are 
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Fig. 2. HRD of the Lupus association. Red squares and blue circles denote, respectively, the WTTSs and CTTSs in our sample. The green 
triangles represent those stars with transition disks (see Sect. 4). The error bars indicate the Icr uncertainties on the stellar luminosities and effective 
temperatures. The black solid lines and dashed lines indicate, respectively, the evolutionary tracks and isochrones computed by Siess et al. (2000) 
with Y = 0.277 and Z = 0.02. The green solid line denotes the zero-age main sequence (ZAMS) computed with that model. The dash-dotted lines 
show the disk lifetime Tj for the TTSs in the Lupus (this work) and Taurus-Auriga (Bertout et al. 2007) associations. 


Table 3. Mean masses and ages of the TTSs in the Lupus association derived from the D’Antona & Mazzitelli (1997) and Baraffe et al. (1998) 
pre-main sequence stars evolutionary models. 



D’Antona & Mazzitelli (1997) 


Baraffe et al. (1998) 


Sample 

Mean 

M/Mq 

Median 

MIMq 

Mean 

logt 

(t in Myr) 

Median 

logf 

(t in Myr) 

Mean 

M/Mq 

Median 

MIMq 

Mean 

logf 

(t in Myr) 

Median 

logf 

(t in Myr) 

CTTS 

0.40 ± 0.02 

0.38 

6.15 ±0.21 

6.10 

0.63 ± 0.04 

0.60 

6.38 ± 0.20 

6.35 

WTTS (on-cloud) 

0.61 ± 0.06 

0.47 

6.79 ± 0.59 

6.42 

0.72 ± 0.06 

0.73 

7.05 ± 0.45 

6.80 

WTTS (off-cloud) 

1.00 ±0.07 

1.00 

6.93 ± 0.35 

6.79 

0.99 ± 0.05 

0.98 

7.17 ±0.28 

7.09 

WTTS (all stars) 

0.80 ± 0.05 

0.80 

6.86 ± 0.46 

6.63 

0.85 ± 0.04 

0.86 

7.11 ±0.36 

7.04 
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Table 4. Results of the Kolmogorov-Smirnov (KS) and Wilcoxon rank sum (WRS) statistical tests applied to the age distributions of the TTSs in 
the Lupus association computed with dilferent sets of isochrones. We provide the p-value given by each test. 


Siess et al. (2000) D’Antona & Mazzitelli (1997) Baraffe et al. (1998) 


Sample 1 

Sample 2 

Pks 

PWRS 

Pks 

PWRS 

Pks 

PWRS 

CTTS 

WTTS (on-cloud) 

2.13 X 10-^ 

1.31 X 10-^ 

3.06 X 10-® 

3.52 X 10-® 

1.65 X 10-^ 

2.71 X 10-® 

CTTS 

WTTS (off-cloud) 

1.01 X 10-'® 

2.01 X 10-® 

5.66 X 10-® 

4.27 X 10-^ 

3.27 X 10-® 

1.24 X 10-® 

CTTS 

WTTS (all stars) 

4.19 X 10-* 

1.57 X 10-® 

3.95 X 10-*' 

2.31 X 10-® 

1.66 X 10-® 

4.25 X 10-® 


Table 5. Results of the Kolmogorov-Smirnov (KS) and Wilcoxon rank sum (WRS) statistical tests applied to the mass distributions of the TTSs 
in the Lupus association computed with different evolutionary tracks. We provide the p-value given by each test. 




Siess et al. (2000) 

D’Antona & Mazzitelli (1997) 

BarafFe et al. (1998) 

Sample 1 

Sample 2 

Pks 

PWRS 

Pks 

PWRS 

Pks 

PWRS 

CTTS 

WTTS (on-cloud) 

3.84 X lO-*' 

2.13 X 10-' 

1.95 X 10-^ 

7.72 X 10-^ 

1.82 X 10"' 

2.95 X 10-' 

CTTS 

WTTS (off-cloud) 

1.42 X 10-^ 

1.10 X 10-® 

9.25 X 10-'® 

4.79 X 10-" 

9.88 X 10-® 

1.62 X 10-® 

CTTS 

WTTS (all stars) 

6.16 X 10-® 

2.39 X 10-'' 

1.97 X 10-® 

4.20 X 10-® 

5.45 X lO-"* 

2.28 X 10-" 


marginally above the adopted significance level in a few cases 
(see Table 5) and remain inconclusive. 

Furthermore, we performed Monte Carlo simulations to as¬ 
sess the robustness of this finding and the errors that affect our 
mass and age estimates. To do so, we constructed a total of 1000 
synthetic samples of the Lupus association by varying, within 
their error bars, the effective temperatures and the remaining ob¬ 
servables (parallaxes, magnitudes, and visual extinction) used 
to compute the stellar luminosities of each star (as described in 
Sect. 2), assuming that these variables are normally distributed. 
Then, we calculated the mass and age of each star from the Siess 
et al. (2000) models and performed a Kolmogorov-Smirnov and 
a Wilcoxon rank sum test on the mass and age distributions of the 
CTTSs and WTTSs for each synthetic dataset. The result of this 
analysis is the mass and age distributions for all synthetic realiza¬ 
tions of the Lupus comoving stars (see Fig. 3), and the histogram 
of probabilities that result from the statistical tests (see Fig. 4). 
It is apparent from Fig. 3 that the on-cloud and off-cloud WTTSs 
define different age distributions despite exhibiting similar kine¬ 
matic properties (as discussed in Paper I). The mass distribution 
of both TTS subclasses shows that there is a lack of CTTSs at 
higher masses, while the WTTSs display a wide range of masses 
(see also Wichmann et al. 1997a). 

To investigate whether the different mass and age distribu¬ 
tions of the TTSs in the Lupus association could be an artefact 
of the procedure used in Sect. 2 to compute the stellar luminosi¬ 
ties, we utilized a different set of tables to convert the magnitudes 
and spectral types to luminosities and effective temperatures. We 
built a control sample by using the bolometric corrections and 
temperatures tabulated by Pecaut & Mamajek (2013) to recom¬ 
pute the stellar luminosities. We then repeated the Monte Carlo 
simulations with our control sample as described above. The re¬ 
sults of the KS- and Wilcoxon statistics applied to the control 
sample are also shown in Fig. 4. These results confirm the dif¬ 
ferent age distributions of the CTTSs and WTTSs in the Lupus 
association for our data and control samples. 

To summarize, the statistical investigations reported above of 
the mass and age distributions of CTTSs and WTTSs contained 
in the Lupus association, using different conversions between 
observed data and stellar parameters, as well as the different evo¬ 
lutionary models, make us confident that these two subgroups 
are indeed different as far as ages are concerned. A similar con¬ 
clusion was reached for the TTS population in the Taurus-Auriga 
association by Bertout et al. (2007). 


4. T Tauri disk iifetimes in the Lupus association 

In their study of the Taurus-Auriga SFR, Bertout et al. (2007) 
obtained a relationship between the disk lifetime and the mass 
of its parent star. We performed a similar study for the Lupus 
association, with the aim of investigating how the different envi¬ 
ronment of the Lupus SFR affects the evolution of protoplanetary 
disks. We thus assumed that evolution from CTTSs to WTTSs is 
caused by the accretion of their circumstellar disks. Following 
Bertout et al. (2007), the disk mass of a TTS is written as 

Md = a(MJM^-f, (2) 

while the mass accretion rate in units of M^lyr is given by 

Mace = , (3) 

where a, jS, and y are unknowns to be determined in our analy¬ 
sis. This last relationship agrees with previous studies (see, e.g., 
Muzerolle et al. 2003, 2005), but more recent work (Alcala et al. 
2014) indicates that the exponent in Eq. 3 might be somewhat 
lower for low-mass CTTSs. We will come back to this point be¬ 
low, but for the time being we use Eq. 3 to compare our findings 
to those of Bertout et al. (2007). The disk lifetime = M^IMacc 
can then be written as 

logTrf = log(a/y) + (J3-2.1) log(M*/Mo). (4) 

We varied these parameters in the range of 5 < log(Q;/ 7 ) < 8 
and -4 < jS < 4 in steps of 0.01 to compute for each star 
with the stellar masses given in Table 1. Then, we compared 
this derived disk lifetime Td for each star with its isochronal age 
f*, also listed in Table 1. The model developed by Bertout et al. 
(2007) predicts that the star is a CTTS whenever f* < Trf(M*), or 
a WTTS if f* > t^(M*). Einally, we compared the mass distri¬ 
butions of the modeled CTTSs (and WTTSs) with the observed 
CTTSs (and WTTSs) in our data sample and ran a KS-test to 
find the best match between these populations for a given com¬ 
bination of a, jS, and y. We adopted a probability threshold of 
99.95% (i.e., p^s > 0.9995) as in Bertout et al. (2007). Doing 
so, we find log(a/ 7 ) = 6.53 + 0.05 and fi - 2.65 + 0.10. Our 
analysis is consistent with an average disk lifetime given by 

Trf ^ 3 X 10® (M*/Mo)“ ®®<^° '°)yr- (5) 

The position of the thus-derived disk lifetime of Lupus TTSs is 
plotted in the HRD of Eig. 2, together with the result obtained 
by Bertout et al. (2007) for the Taurus-Auriga SER. 
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Fig. 3. Mass and age distributions of the CTTSs and WTTSs of the Lupus association obtained after 1000 Monte Carlo simulations. The green 
solid line indicates the kernel density estimator. 


To assess the significance of this finding, we performed 
Monte Carlo simulations. We applied the above strategy to the 
1000 synthetic realizations of the Lupus association, as well as 
to the control sample constructed in Sect. 3. The results, shown 
in Table 6, demonstrate that the simulations are consistent and 
confirm the values mentioned above. Simulations of the control 
sample returned lower values for the disk parameters, but these 
are statistically compatible with our result (within 3cr). Thus, 
our results remain valid if we use a different set of photomet¬ 
ric tables and spectral type to effective temperature calibration 
for computing the stellar luminosities and the resulting mass and 
age estimates. The evolution of TTSs of the Lupus association is 
therefore consistent with the evolutionary scenario proposed by 
Bertout et al. (2007) for the Taurus-Auriga T association. 

As mentioned above, Alcala et al. (2014) have recently found 
that Macc in a sample of Lupus low-mass accreting 

stars. If we use this revised accretion rate estimate in Eqs. 3 and 
4 and run our calculations as described above, the previous value 


of P - 2.65 + 0.10 shifts to p - 2.35 + 0.10, and the resulting 
exponent on the stellar mass in Eq. 5 is not affected. Similarly, 
if we use the mass accretion rate of Macc derived by 

Venuti et al. (2014) for a sample of CTTSs in the NGC 2264 open 
cluster, we find that p - 2.65 + 0.10 shifts to p - 1.95 + 0.10. In 
both cases, the different value obtained for p affects the disk mass 
given by Eq. 2 in our disk model. However, the combination of 
the different exponent p on the disk mass with a different mass 
accretion rate Macc yields the same exponent on the disk lifetime 
in Eq. 5. Thus, our result for the average disk lifetime remains 
the same. 

Recently, Daemgen et al. (2013) have investigated the impact 
of binary companions on the evolution of circumstellar disks and 
concluded that they evolve faster in binary systems. In this con¬ 
text, we decided to remove the binaries (and multiple systems) 
from our sample to inspect this point. Doing so, we recalculated 
the disk lifetime in Lupus and find log(a/y) = 6.56 + 0.03 and 
P - 2.69 + 0.06. Thus, the average disk lifetime excluding bina- 
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Data Sample: Kolmogorov-Smirnov Test 


Data Sample: Wilcoxon Test 




Control Sample: Kolmogorov-Smirnov Test 
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Fig. 4. Histogram of prohahilities (i.e., p-value) obtained from the Kolmogorov-Smirnov {left panels) and Wilcoxon tests {right panels) for 
the age distributions of the CTTSs and WTTSs in our data sample {upper panels) and control sample {lower panels) after 1000 Monte Carlo 
realizations. The red dashed line indicates the significance level of cr = 0.05. 


ries from the sample is ^ 4 x 10^ This 

result implies that the average disk lifetime of the association is 
likely to be shorter when binaries (and multiple systems) are in¬ 
cluded in our analysis. However, it is important to mention that 
both results are still statistically compatible within their error 
bars. For a more direct comparison with the results obtained by 
Bertout et al. (2007) in the Taurus-Auriga region (see Sect. 5), 
we prefer our first solution, because the binary stars were not 
excluded in their study. 


Table 6. Mean and median values of the disk parameters {a, fj, and y) 
obtained via Monte Carlo simulations after 1000 realizations of the data 
and control samples. 



log(a/r) 

ys 



Mean Median 

Mean 

Median 

Data sample 

6.52 + 0.05 6.52 

2.51 + 0.10 

2.52 

Control sample 

6.35 ± 0.05 6.36 

2.21 + 0.09 

2.22 


Despite the success of this empirical disk model in explain¬ 
ing the different age distributions of the CTTSs and WTTSs in 
the Lupus association, we stress that the conclusions are only 
meaningful in a statistical sense and that our model may lead to 
a false classification in a few cases. Three WTTSs (RXJ1559.9- 
3750, RXJ1609.4-3850, and RXJ1612.0-3840) were misidenti- 
hed as CTTSs in our analysis because of an inconsistency be¬ 
tween the observed isochronal age and the computed lifetime 
of their disks at the Icr level. These misidentified WTTSs are 
among the youngest stars in our sample. 

Another interesting point in our study is the existence of 
six stars (Sz91, Sz95, Sz96, Szll2, RXJ1609.9-3923, and 
RXJ1608.5-3847) in our sample that were first identihed as 
CTTSs and later reclassified as transition disk objects (see 
Romero et al. 2012; Tsukagoshi et al. 2014). We have labeled 
them with the mention "CTTS/TD" in Table 1 to distinguish 
them from the remaining CTTSs. Romero et al. (2012) claim 
that RXJ1609.9-3923, the youngest transition disk object in the 
sample, is a “triple system with tight binary components consis¬ 
tent with two equally bright objects". Later in the text, they state 
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that the source has an uncertain classification because its spectral 
energy distribution resembles a CTTS of M spectral type. Taking 
only the position of RXJ1609.9-3923 into account in the HRD, 
our analysis suggests that this source is more likely a CTTS. 
The remaining transition disk objects in the sample intermingle 
among CTTSs and WTTSs in the HRD, and their positions are 
fully consistent with the disk lifetime locus that roughly sepa¬ 
rates both species (see Fig. 2). 

5. Comparison with the Taurus-Auriga association 

Both the Lupus and Taurus-Auriga associations discussed in this 
paper have been identified by using the convergent point search 
method applied to proper motion data in our previous studies 
(Bertout & Genova 2006; Galli et al. 2013) and the physical 
properties of the comoving stars were computed in the same 
completeness limits as described in Sect. 2. In this context, 
it is interesting to note that the Lupus association seems to be 
older than the Taurus-Auriga association. It is also apparent 
(see Fig. 2) that the distribution of CTTSs in Lupus is shifted 
toward lower masses and late spectral types, while in Taurus- 
Auriga they display a wider range of masses and effective tem¬ 
peratures (see, e.g.. Fig 1 of Bertout et al. 2007). These first 
results had already been suggested in previous studies by us¬ 
ing an average distance to the Lupus SFR to compute the stellar 
ages (and masses) for all stars and comparing different samples 
of TTSs in these SFRs (Hughes et al. 1994; Wichmann et al. 
1997a). The individual parallaxes from Paper I allowed us to 
unambiguously confirm this scenario by comparing two samples 
of TTSs that have been investigated with the same selection cri¬ 
teria. 

Bertout et al. (2007) derived an average disk lifetime of ^ 
4x10^ (M*/MO)°'^^ yr for the TTSs in the Taurus-Auriga associ¬ 
ation that contrasts with the result of ^ 3x 10^ (M*/Mo)°'^^ yr 
obtained in this paper for the Lupus association. Our heuristic 
disk model thus predicts different disk lifetimes for the two as¬ 
sociations. First, the exponents defining the mass dependence 
are different, but the two values are still compatible with each 
other if one considers the uncertainties involved {[5 — 2.65+0.10 
in Lupus and fi — 2.85 + 0.21 in Taurus-Auriga). Second, re¬ 
gardless of the exact p value, we note that the standard deviation 
values of log(a/y) found from the Monte Carlo simulations of 
Lupus and Taurus-Auriga comoving stars lead to disjointed age 
distributions at the Icr significance level for the IMg stars of 
both regions (2.95 - 3.72 Myr for Lupus and 3.80 - 4.79 Myr 
for Taurus-Auriga). Although this result will need confirmation 
with more significant samples statistically, it gives a first indica¬ 
tion that disk lifetimes can be different in different star-forming 
regions. 

This is not the only indication that disk lifetimes are likely to 
be different in Taurus-Auriga and Lupus. For example, a com¬ 
parative analysis of the HRD of each association reveals that the 
overlap between CTTSs and WTTSs is more important in Lu¬ 
pus than in Taurus-Auriga (see Fig. 2 of this paper and Fig. 1 
of Bertout et al. 2007). According to the results for the Taurus- 
Auriga association obtained by Bertout et al. (2007), it is only 
the evolutionary status of the circumstellar disks that allows us 
to distinguish between both subgroups. In this context, the exist¬ 
ing overlap between both species in the HRD could be indicative 
of a shorter disk lifetime for the TTS population in the Lupus as¬ 
sociation as illustrated in Fig. 2. 

Another fact that seems to support the different disk life¬ 
times in each association is the age of the transition disk ob¬ 
jects that show signs of disk dissipation and are presumably in 


an intermediate state between CTTSs and WTTSs with little or 
no detectable near-IR excess emission and significant far-IR ex¬ 
cess (Calvet et al. 2005; Romero et al. 2012). When excluding 
RXJ1609.9-3923 for the reasons discussed in Sect. 4, the aver¬ 
age age of the remaining transition disks included in our sam¬ 
ple (Sz91, Sz95, Sz96 , Szll2, and RXJ1608.5-3847) is only 
1.9 + 0.3 Myr. On the other hand, the average age of the transi¬ 
tion disks (DM Tau, GM Aur, LkCa 15, and DI Tau) in Taurus- 
Auriga is 5.3 + 1.4 Myr (see Table 4 of Bertout et al. 2007), thus 
significantly older than the observed value for Lupus. 

Martin et al. (1998) investigated the ratio of WTTSs over 
CTTSs in different regions of the p Ophiuchi molecular cloud 
complex, and they argue that the high WTTS/CTTS ratio ob¬ 
served in some fields could be the result of the strong UV radia¬ 
tion, stellar winds, or supernova explosions from nearby massive 
stars that act to shorten the lifetime of circumstellar disks (see 
also Walter et al. 1994). A similar argument could apply to the 
Lupus association, where the WTTSs greatly exceed the number 
of CTTSs as compared to the Taurus-Auriga association. The 
Lupus cloud complex is located in an environment disturbed by 
the existence of hot stars from the Scorpius-Centaurus OB asso¬ 
ciation and their dynamical effects, while the low-mass SFR of 
Taurus-Auriga is more isolated and shows no (or weak) interac¬ 
tion with hot stars. In this context, the large number of WTTSs 
in Lupus could therefore be caused by the early disappearance 
of the circumstellar disks in CTTSs. 

6. Conclusions 

The newly derived individual parallaxes from our previous kine¬ 
matic study of the Lupus association have been used in this paper 
to determine the photospheric luminosities and refine the masses 
and ages of the TTS population in this SFR. We investigated the 
mass and age distributions of CTTSs and WTTSs in the Lupus 
association and demonstrated that the probability of both sub¬ 
groups being drawn from the same age distribution is very low. 
We also confirmed the different properties of on-cloud and off- 
cloud WTTSs of the Lupus association using our individual par¬ 
allaxes, which is one result that had been anticipated by Wich¬ 
mann et al. (1997a) by using an average distance estimate of the 
Lupus region for all stars. We concluded that the different mass 
and age distributions in the Lupus association can be explained 
by assuming that a CTTS evolves into a WTTS when the disk is 
accreted by the central star. Such a result had already been pro¬ 
posed by Bertout et al. (2007) for the Taurus-Auriga association, 
and it was now confirmed for a different SFR. 

Based on an empirical disk evolution model, we derived the 
average lifetime of a circumstellar disk in the Lupus associa¬ 
tion in terms of the stellar mass and compared our result with a 
similar study conducted for the Taurus-Auriga association. We 
argued that specific properties of the Lupus association, includ¬ 
ing the more evident overlap of CTTSs and WTTSs in the HRD, 
the age of transition disk objects, and the ratio of WTTS/CTTS, 
could be the result of an earlier disappearance of the circumstel¬ 
lar disks in CTTSs implying that the timescale of disk dissipation 
in the Lupus association is shorter (than in Taurus-Auriga). 

This study represents an important step toward a better un¬ 
derstanding of the early stages of star formation and disk evo¬ 
lution in different environments. It will be extended to other 
nearby SFRs in the future to allow for a more detailed compara¬ 
tive analysis of the evolution of TTSs in different associations. 
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Table 1. Physical properties of the Lupus comoving stars. 


Star 

ST 

Teff 

EW(Ha) 

Ay 

L/Lo 

Filter 

M/Mo 

logt 

T Tauri 

logT 

Ref. 



(K) 

(A) 

(mag) 




(t in Myr) 

Subclass 

(r in Myr) 


RXJ1508.8-3715 

K5 

4350 

1.02 

0.25 ±0.17 

1 24+1-f^O 

43073:5? 

0257i:i 

0 36+^-^^ 
-Q16 

0 31+0.45 

0 547-^4 

V 

0.80 ±0.15 

7.59 ± 0.26 

WTTS 

6.48 ±0.14 

1,2 

RXJ1518.4-3738 

KO 

5250 

1.91 

0.32 ±0.10 

Ic 

0.92 ±0.18 

7.52 ±0.19 

WTTS 

6.51 ±0.15 

1.2,3 

RXJ1524.5-3652 

KO 

5250 

1.35 

0.44 ± 0.07 

Ic 

1.17±0.27 

7.29 ± 0.32 

WTTS 

6.57 ±0.17 

1,2,3 

RXJ1525.0-3604 

K1 

5080 

0.94 

0.53 ±0.18 

Ic 

1.68 ±0.45 

6.77 ± 0.39 

WTTS 

6.65 ± 0.20 

1,2,3 

RXJ1525.5-3613 

K2 

4900 

0.59 

0.51 ±0.21 

Ic 

1.29 ±0.34 

7.06 ± 0.42 

WTTS 

6.59 ± 0.20 

1,2,4 

RXJ1531.3-3329 

G9 

5410 

1.68 

1.54 ±0.04 

Ic 

1.73 ±0.54 

6.83 ± 0.34 

WTTS 

6.66 ± 0.23 

1,2,3 

RXJ1534.6-4003 

K5 

4350 

0.78 

0.33 ±0.16 

Ic 

0.82 ±0.16 

7.50 ±0.36 

WTTS 

6.48 ±0.15 

1.2,3 

RXJ1540.7-3756 

K6 

4205 

0.93 

0.00: 

Ic 

0.89 ±0.15 

7.05 ± 0.58 

WTTS 

6.50 ±0.13 

1,2,5 

RXJ1544.5-3521 

K5 

4350 

0.51 

0.59 ±0.19 

V 

0.90 ± 0.25 

7.36 ±0.39 

WTTS 

6.50 ±0.21 

1,2 

Sz73 

MO 

3850 

97.20 

2.71 ±0.11 

Ic 

0.57 ± 0.08 

6.34 ±0.58 

CTTS 

6.40 ±0.11 

2,6 

GQLup* 

K7-M0 

3955 

38.60 

0.00: 

2 Cl +3.96 

1 927i:3? 

0397i:3? 

023^8:15 

-0-Q9 

0 72+2.20 

-1 12 

Ic 

0.73 ±0.13 

5.82 ±0.35 

CTTS 

6.45 ±0.13 

2,6,7 

RXJ1549.9-3629 

K2 

4900 

1.06 

0.28 ±0.18 

Ic 

1.55 ±0.55 

6.80 ± 0.52 

WTTS 

6.63 ± 0.26 

1,2,4 

RXJ1552.3-3819 

K7 

4060 

0.68 

0.00: 

Ic 

0.80 ±0.12 

6.87 ± 0.54 

WTTS 

6.48 ±0.12 

1,2,4 

RXJ 1605.7-3905 

KO 

5250 

0.30 

0.15 ±0.07 

Ic 

1.25 ±0.46 

7.20 ± 0.46 

WTTS 

6.58 ± 0.27 

2,3 

F304 

K6 

4205 

1.26 

0.00: 

Ic 

0.79 ±0.15 

7.42 ± 0.36 

WTTS 

6.47 ±0.14 

1,2,4 

RXJ1608.5-3847 

M2 

3580 

5.98 

1.00 ±0.07 

Ic 

0.40 ± 0.05 

6.04 ± 0.23 

CTTS/TD 

6.31 ±0.11 

1,2,4,8 

RXJ1610.0-4016 

K2 

4900 

0.88 

0.52 ± 0.28 

Ic 

1.59 ±0.57 

6.76 ±0.51 

WTTS 

6.64 ± 0.26 

1,2,3 

SZ121 

M3 

3470 

6.50 

0.98 ±0.18 

0.52*>«3f 

Ic 

0.31 ±0.08 

6.17 ±0.31 

WTTS 

6.25 ± 0.20 

2,6 

RXJ1613.0-4004 

K7 

4060 

0.74 

1.34 ±0.04 

0-26131! 

Ic 

0.80 ±0.13 

7.15 ±0.71 

WTTS 

6.48 ±0.13 

1,2,4 

RXJ1511.0-3252AB 

K6 

4205 

0.62 

0.00: 

n oo+8:^t 
0.83 „ 

0 39+0-2“ 
-0.13 

0.31+°'|’ 

Ic 

0.94 ±0.12 

6.53 ± 0.43 

WTTS 

6.51 ±0.11 

1,2,4 

RXJ1511.6-3550 

K5 

4350 

0.58 

0.00: 

Ic 

0.96 ±0.15 

7.23 ± 0.40 

WTTS 

6.52 ±0.12 

1.2,4 

RXJ1512.6-3417 

K5 

4350 

_ 

0.56 ±0.19 

V 

0.90 ±0.19 

7.36 ± 0.44 

WTTS 

6.50 ±0.16 

2,9 

HD135127 

F5 

6440 

4.13 

0.10 ±0.09 

3 

4o7l:5§ 

05>5 

aUl 

Ic 

1.41 ±0.08 

7.30 ± 0.20 

WTTS 

6.61 ± 0.06 

1,2,4 

RXJ1515.7-3332 

KO 

5250 

1.68 

0.37 ±0.18 

Ic 

1.22 ±0.30 

7.22 ± 0.36 

WTTS 

6.58 ±0.19 

1,2,3 

RXJ1525.6-3537 

K6 

4205 

0.14 

0.24 ±0.15 

Ic 

0.95 ±0.11 

6.88 ± 0.54 

WTTS 

6.52 ±0.10 

1,2,4 

RXJ1527.3-3603 

K7-M0 

3955 

0.12 

0.50 ±0.21 

V 

0.70 ±0.11 

6.98 ± 0.60 

WTTS 

6.45 ±0.12 

1,2 

RXJ1529.3-3737 

M3 

3470 

0.24 

0.00: 

0.36!; 

05>:i 

Ic 

0.35 ± 0.06 

6.24 ±0.19 

WTTS 

6.28 ±0.14 

1,2,4 

RXJ1529.7-3628 

K2 

4900 

0.89 

0.46 ± 0.20 

V 

1.54 ±0.34 

6.81 ±0.41 

WTTS 

6.63 ±0.17 

1,2,3 

Sz65* 

MO 

3850 

19.40 

0.29 ±0.11 

Ic 

0.56 ± 0.08 

6.00 ± 0.36 

CTTS 

6.39 ±0.12 

2,6,10 

RXJ1539.7-3450 

K4 

4590 

0.16 

0.50 ±0.21 

V 

1.04 ±0.25 

7.24 ± 0.49 

WTTS 

6.54 ±0.18 

1,2 

RXJ1540.3-3426A 

M3.5 

3420 

3.03 

0.00: 


Ic 

0.31 ±0.05 

6.48 ± 0.27 

WTTS 

6.25 ±0.14 

1,2,11 

SSTc2dJ154148.3-350145 

M2.7 

3500 


0.00: 


Ic 

0.36 ± 0.07 

6.45 ± 0.30 

WTTS 

6.28 ±0.15 

2,11,12 

RXJ1542.0-3601 

K7 

4060 

0.32 

0.23 ±0.18 

0 20+8-38 

1 6278:ie 

° -0-83 
0.29+°'® 

V 

0.76 ±0.12 

7.34 ± 0.44 

WTTS 

6.46 ±0.12 

1,2 

RXJ 1544.0-3311* 

KO 

5250 

0.87 

0.65 ±0.19 

V 

1.33 ±0.36 

7.11 ±0.37 

WTTS 

6.60 ± 0.20 

1,2,3 

RXJ1546.6-3618 

K1 

5080 

0.67 

0.58 ±0.33 

Ic 

1.38 ±0.46 

7.03 ± 0.45 

WTTS 

6.61 ±0.25 

1,2,3 

RXJ 1546.7-3459 

MO 

3850 

1.53 

0.00: 

J 

0.59 ± 0.09 

6.74 ± 0.49 

WTTS 

6.40 ±0.12 

1,2 

RXJ1547.1-3540 

K5 

4350 

0.55 

0.20 ±0.18 

0 . 24 + 0:24 

V 

0.80 ±0.10 

7.61 ±0.31 

WTTS 

6.48 ±0.10 

1,2 

RXJ1547.6-4018 

K1 

5080 

1.60 

0.35 ± 0.07 

i.osJI 

-Q.IQ 

1 27+2.67 

034^8:31 

Ic 

1.17 ±0.28 

7.24 ± 0.33 

WTTS 

6.57 ±0.18 

1,2,3 

HMLup 

M3 

3470 

115.30 

0.00: 

Ic 

0.34 ± 0.08 

6.37 ± 0.27 

CTTS 

6.27 ±0.17 

2,6 

HNLup* 

M1.5 

3650 

49.60 

2.14±0.11 

Ic 

0.42 ± 0.09 

5.88 ± 0.25 

CTTS 

6.32 ±0.17 

2,6,10 

RXJ1548.1-3452 

M2.5 

3525 

2.75 

0.00: 

Ic 

0.37 ± 0.05 

6.29 ± 0.24 

WTTS 

6.29 ±0.12 

1,2,4 

RXJ1548.9-3513 

K6 

4205 

0.16 

0.51 ±0.18 

0.15!°‘5 

V 

0.70 ±0.11 

7.79 ± 0.44 

WTTS 

6.44 ±0.12 

1,2 

Sz76 

Ml 

3720 

10.30 

0.38 ±0.18 

019+8:32 

Ic 

0.47 ± 0.05 

6.80 ± 0.55 

WTTS 

6.35 ±0.10 

2,6,8 

HD141277* 

K1 

5080 

1.52 

0.49 ± 0.34 

1 

062^8:35 

-0.12 

0.47™2 

Ic 

1.57 ±0.54 

6.87 ± 0.46 

WTTS 

6.64 ± 0.25 

1,2,3 

RXJ1550.7-3828 

K7 

4060 

0.88 

0.86 ± 0.30 

Ic 

0.75 ±0.12 

6.20 ± 0.47 

WTTS 

6.46 ±0.13 

1,2,4 

Sz77* 

MO 

3850 

12.40 

0.00: 

Ic 

0.57 ± 0.08 

6.27 ± 0.38 

CTTS 

6.39 ±0.12 

2,6,10 

RXJ1555.4-3338 

K5 

4350 

— 

0.00: 

Ic 

0.91 ±0.15 

7.35 ± 0.41 

WTTS 

6.51 ±0.13 

1,2,4 

RXJ 1556.0-3655 

Ml 

3720 

73.60 

0.14 ±0.07 

Ic 

0.48 ± 0.05 

6.28 ± 0.37 

CTTS 

6.35 ±0.10 

1,2,4 

Sz82 * 

MO 

3850 

39.00 

0.00: 

0.93!8:^4 

Ic 

0.56 ± 0.08 

6.06 ± 0.35 

CTTS 

6.39 ±0.12 

2,6,8,10 

Szl26* 

K7-M0 

3955 


0.60 ±0.13 

0 43+^-^^ 
0.36™^ 

Ic 

0.68 ±0.12 

6.66 ± 0.54 

CTTS 

6.44 ±0.14 

2,6,10,13 

Szl28* 

M1.5 

3650 

9.90 

1.86 ±0.13 

Ic 

0.43 ±0.10 

6.36 ± 0.41 

CTTS 

6.33 ±0.17 

2,6,14 

RXJ1558.9-3646 

M1.5 

3650 

2.81 

0.00: 

0 46+^-36 

Ic 

0.43 ±0.10 

6.24 ±0.31 

WTTS 

6.33 ±0.17 

1,2,4 

CD-3610569 

K3 

4730 

0.71 

0.58 ± 0.08 

0 69+^-^'^ 

015^8:1^ 

-Q.Q5 

0.55!g 

0.88!|| 

1 09+0.85 

1 2Q+0.94 

-Q.0§ 

0.57!«; 

0.80!j; 

'^•■^^- 0.15 

0 93+1-04 

0 40+^-^^ 

^•^^-0 17 

Ic 

1.10 ±0.27 

7.24 ± 0.39 

WTTS 

6.55 ±0.18 

1,2,3 

RXJ1559.9-3750 

M0.5 

3785 

0.51 

0.32 ± 0.28 

Ic 

0.51 ±0.10 

5.98 ± 0.47 

WTTS 

6.37 ±0.15 

1,2,4 

SSTc2dJ 160000.6-422158 

M3.7 

3400 

— 

0.00: 

Ic 

0.30 ± 0.05 

6.57 ±0.31 

WTTS 

6.24 ±0.13 

2,11,12,15 

Szl31 

M2 

3580 

31.90 

2.48 ± 0.24 

Ic 

0.41 ± 0.06 

6.12 ±0.77 

CTTS 

6.31 ±0.12 

2,6,15 

RXJ 1601.9-3613 

K3 

4730 

0.45 

0.59 ±0.21 

Ic 

1.20 ±0.29 

7.10 ±0.48 

WTTS 

6.57 ±0.18 

1,2,4 

EXLup 

MO 

3850 

43.30 

0.00: 

Ic 

0.56 ± 0.08 

5.99 ±0.31 

CTTS 

6.39 ±0.12 

2,6 

RXJ1603.8-3938* 

K3 

4730 

1.13 

0.43 ± 0.08 

Ic 

1.01 ±0.25 

7.35 ± 0.39 

WTTS 

6.53 ±0.19 

1,2,3,16 

HD143978 

G2 

5860 

3.60 

0.00: 

Ic 

1.11 ±0.07 

7.88 ± 0.48 

WTTS 

6.56 ± 0.07 

1,2,4 

RXJ1605.5-3837 

Ml 

3720 

1.28 

0.00: 

Ic 

0.50 ± 0.05 

6.74 ± 0.42 

WTTS 

6.36 ± 0.09 

1,2,4 

HOLup* 

Ml 

3720 

219.80 

0.00: 

Ic 

0.48 ± 0.05 

6.19 ±0.31 

CTTS 

6.35 ±0.10 

2,6,10 

Sz90* 

K7-M0 

3955 

28.50 

2.12 ±0.40 

Ic 

0.66 ±0.13 

6.26 ± 0.65 

CTTS 

6.43 ±0.15 

2,6,10 

Sz91* 

M0.5 

3785 

95.90 

1.93 ±0.11 

Ic 

0.52 ±0.10 

6.51 ±0.44 

CTTS/TD 

6.38 ±0.15 

2,6,10 

RXJ1607.2-3839 

K7 

4060 

1.40 

0.29 ± 0.20 

Ic 

0.75 ±0.11 

6.29 ± 0.47 

WTTS 

6.46 ±0.12 

1,2,4 

Sz95 * 

M1.5 

3650 

10.20 

1.29 ±0.11 

Ic 

0.43 ±0.10 

6.30 ± 0.37 

CTTS/TD 

6.33 ±0.17 

2,6,15 
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Table 1. continued. 


Star 

ST 

Teff 

EW{Ha) 

Ay 

L/Lo 

Filter 

M/Mo 

logt 

T Tauri 

logT 

Ref. 



(K) 

(A) 

(mag) 




(t in Myr) 

Subclass 

(r in Myr) 


RXJ1608.0-3857 

MO 

3850 

1.14 

0.24 ± 0.22 

092+™ 

05v8:^^ 

0-39!oi! 

Ic 

0.56 ± 0.08 

6.06 ± 0.42 

WTTS 

6.39 ±0.12 

1,2,4 

Sz96* 

M1.5 

3650 

11.00 

0.43 ±0.11 

Ic 

0.44 ±0.10 

6.17 ±0.33 

CTTS/TD 

6.33 ±0.17 

2 ,6,10 

RXJ1608.3-3843 

K7 

4060 

0.11 

0 .00: 

Ic 

0.80 ±0.12 

6.87 ± 0.44 

WTTS 

6.48 ±0.12 

1,2,4 

Sz97 

M3 

3470 

58.20 

0.95 ±0.11 

0.28+11 

Ic 

0.34 ± 0.08 

6.34 ± 0.28 

CTTS 

6.28 ±0.17 

2,6 

RXJ1608.6-3922 

K6 

4205 

7.42 

1.52 ±0.07 

0 .88f™ 

0.51+°” 

Ic 

0.93 ±0.13 

6.49 ± 0.47 

CTTS 

6.51 ±0.11 

1 ,2,4,8 

SSTc2dJ160853.2-391440 

K7.8 

3900 

_ 

4.60 ± 0.43 

Ic 

0.61 ±0.10 

6.45 ± 0.70 

CTTS 

6.41 ±0.13 

2,11,12,15 

RXJ1608.9-3905 

K2 

4900 

0.84 

0.38 ±0.18 

oi+Jf 

Ic 

1.86 ±0.34 

6.47 ± 0.40 

WTTS 

6.68 ±0.14 

1,2,4 

RXJ1608.9-3945 

M3.5 

3420 

2.69 

0.16 ±0.14 

Ic 

0.32 ± 0.06 

6.36 ± 0.28 

WTTS 

6.26 ±0.15 

1,2,4 

Szlll 

M1.5 

3650 

145.20 

0 .00: 

0.34f|; 

-QQ6 
0-25 onQ 

Ic 

0.43 ±0.10 

6.38 ±0.33 

CTTS 

6.33 ±0.17 

2,6 

Szll2 

V908SCO 

M4 

M4 

3370 

3370 

46.70 

33.00 

0.81 ±0.11 
0 .00: 

Ic 

Ic 

0.30 ± 0.05 
0.30 ± 0.04 

6.26 ± 0.29 
6.18 ±0.55 

CTTS/TD 

CTTS 

6.24 ±0.14 
6.24 ±0.13 

2,6 

2,6 

Szll5 

M4 

3370 

7.10 

0 .00: 

Ic 

0.29 ± 0.05 

6.55 ± 0.30 

WTTS 

6.23 ±0.15 

2,6,15,19 

Szl34 

Ml 

3720 

90.50 

0 .00: 

Ic 

0.47 ± 0.05 

6.64 ± 0.40 

CTTS 

6.35 ±0.10 

2,6 

RXJ1609.4-3850 

M0.5 

3785 

1.64 

0.13±0.12 

I4i;?:§^ 

030+8:^^ 

-Q-I5 

0 

0 2^5:^^ 

1 +l+i+l-24 

0 978:51 

-Q-18 

0.26+^-|| 

Ic 

0.51 ±0.08 

5.86 ±0.32 

WTTS 

6.37 ±0.13 

1,2,4 

Szll6* 

M1.5 

3650 

3.90 

0.42 ± 0.24 

Ic 

0.43 ±0.10 

6.45 ± 0.43 

WTTS 

6.33 ±0.18 

2,6,15 

Szll7 

M2 

3580 

20.50 

0.43 ±0.11 

Ic 

0.39 ± 0.05 

6.38 ±0.32 

CTTS 

6.31 ±0.12 

2,6 

Szll8 

K6 

4205 

28.30 

4.19±0.11 

Ic 

0.96 ±0.14 

6.73 ± 0.88 

CTTS 

6.52 ±0.11 

2,6 

RXJ1609.9-3923* 

M1.5 

3650 

13.67 

2.33 ± 0.07 

Ic 

0.42 ± 0.09 

5.88 ±0.28 

CTTS/TD 

6.32 ±0.17 

1,2,4 

Szll9 

M4 

3370 

4.50 

0.53 ±0.18 

Ic 

0.30 ± 0.05 

6.18 ±0.61 

WTTS 

6.25 ±0.13 

2,6 

Szl22 

M2 

3580 

4.80 

0.37 ±0.18 

Ic 

0.39 ± 0.05 

6.45 ± 0.38 

WTTS 

6.31 ±0.12 

2,6 

Szl23* 

Ml 

3720 

250.60 

0.86 ±0.11 

0 39+o'4i 

Ski 

Ic 

0.48 ± 0.05 

6.38 ± 0.43 

CTTS 

6.35 ±0.10 

2 ,6,10 

RXJ1612.0-3840 

K5 

4350 

0.72 

0.61 ±0.26 

Ic 

1.15 ±0.25 

5.87 ±0.41 

WTTS 

6.56 ±0.16 

1,2,4 

SSTc2dJ161207.6-381324 

M4.5 

3305 

3.60 

0 .00: 

2.i9);+'°° 

Ic 

0.25 ± 0.05 

6.14 ±0.64 

WTTS 

6.20 ±0.15 

2,11,12,18 

SSTc2dJ161243.8-381503 

K4.8 

4400 

_ 

1.27 ±0.12 

Ic 

1.20 ±0.18 

6.18 ±0.45 

CTTS 

6.57 ±0.12 

2,11,12,15 

RXJ1614.4-3808 

K7 

4060 

0.46 

0.23 ±0.17 

0 15+^'*^ 

V 

0.70 ±0.11 

7.54 ±0.35 

WTTS 

6.45 ±0.12 

1,2 

HD147402 

G3 

5830 

— 

0.37 ± 0.28 

Ic 

1.14 ±0.16 

7.40 ±0.16 

WTTS 

6.56 ±0.11 

2,3,8 


Notes. For each star, we provide the spectral type, temperature, Ha equivalent width, visual extinction, luminosity, photometric filter used to 
compute luminosities, the mass and age derived from the Siess et al. (2000) models, disk lifetime computed from Eq. (5), T Tauri subclass, and 
sources of photometric/spectroscopic information. The Ay estimate of 0.00: denotes a non-physical negative value or a non-significant result. 
Stars marked with are binary (or multiple) systems (see Sect. 2). 

References. (1) Krautter et al. (1997); (2) 2MASS (Cutri et al. 2003); (3) Torres et al. (2006); (4) Wichmann et al. (1997a); (5) Sartori et al. (2003); 
(6) Hughes et al. (1994); (7) Neuhiiuser et al. (2008); (8) Wahhaj et al. (2010); (9) Wichmann et al. (1997b); (10) Ghez et al. (1997); (11) Comeron 
et al. (2009); (12) Lopez Marti et al. (2011); (13) Valenti et al. (2003); (14) Brandner et al. (1996); (15) Merin et al. (2008); (16) Guenther et al. 
(2001); (17) Alcala et al. (2014); (18) Comeron et al. (2013); (19) Cieza et al. (2007); 
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